Separation of variables for Gaudin-Calogero systems 



B. Enriquez, B. Feigin and V. Rubtsov 

Abstract. We construct an elliptic analogue of Sklyanin's separation of vari- 
ables for the sl(2) Gaudin system, using an adaptation of Drinfeld's Radon transfor- 
mations. 

Introduction. 

The geometric Langlands conjectures, as formulated by Beilinson and Drinfeld in 
[2], aim at the construction of certain D-modules on the moduli space of G-bundles 
over (punctured) curves (G a reductive group). Positive characteristic versions of 
these conjectures had been solved earlier by Drinfeld in [3], in the case G = GLi. 

In the paper [6] , the P-modules arising from the construction of [2] were studied 
in the special case of a rational curve with marked points, and identified with the 
Gaudin model. Then, in [8], Drinfeld's construction of local systems on the moduli 
space of rank two vector bundles on a curve in positive characteritic ([3]) was adapted 
to the complex situation. The identification of these two constructions amounts to 
Sklyanin's separation of variables ([10]), as it was noticed in [8]. This computation 
is recalled in the first part of this text. 

The question has been raised in [8] to construct a similar separation of variables 
for the Gaudin-Calogero systems, which were computed in [4] and [9], and play a 
similar role in the case of a punctured elliptic curve. This note aims at solving this 
question. In the present case Drinfeld's diagrams for Radon transformation have to 
be slightly modified. 

It is also worth to note that the systems presented here, are the specialization 
at the critical level, of the Knizhnik-Zamolodchikov-Bernard equations on the torus. 
According to the general viewpoint that the diagrams of [3] are related to the Drinfeld- 
Sokolov reduction, the generalisation of the present work to noncritical level should 
relate these equations to the Virasoro correlators on the torus. In [5], the Bethe 
equations were connected with the unitarity property of the KZB equation on the 
torus. It would be desirable to clarify further the connection between these issues. 

We express our thanks to A. Stoyanovsky and A. Varchenko for having discussed 
with us the content of this paper, and to C. Sabbah for consultations on V- modules. 
V.R. was supported by the INTAS grants 93-2494 and 1010-CT93-0023, and by the 
CNRS; he expresses his thanks to these institutions. 
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1. Separation of variables for sl 2 Gaudin systems. 

Let us begin with some reminders on the Gaudin system. Let X = CP 1 , z a be 
marked points on X, a = 1, • • • , N, G be >SL2(C), B C G be the upper triangular sub- 
group. The moduli space M.q{X, z a ) of G-bundles on CP 1 , with parabolic structures 

at z a , is the disjoint union of the J\4q (X, z a ), n > 0, corresponding to the parabolic 

structures on the sheaf 0(noo) © 0(—noo). We then identify M.^(X, z a ) with 

P n \(G/B) N , P = G, and P n = P ^ ,t e C x ,p(z) E C[z],deg(p) < 2n} for 

n > 0, Po acting diagonally and P n acting on the ct-th factor by left translation, after 
the replacement of z bj z a . In what follows we will deal with A4q (X, z a ). 

Let us fix weights, A a , a = 1,---,N. On A4q (X, z a ) lives the bundle £(A a ) 5 
quotient of the bundle ^^ =1 C\ a (£\ a is the line bundle on G/B, corresponding to the 
weight A a ). The natural action of Z{U-2sl2)i oc (the center of the local completion 
of the enveloping algebra of the central extension of sl2(C((z))) at level —2) is by 
differential operators L a , a = 1, • • • , N, which were identified in [6] with the Gaudin 

hamiltonians L a = Y\ n ^ n , / a , I a being an orthogonal basis of sl 2 (C). 

Following the conjectures of [2], the ©-modules on M.q(X, z a ) (twisted by £(a„)) 
defined by L a — ^ a should satisfy the Hecke eigenvalue property. In [3], a construction 
of such modules was given in the case of a curve of genus > without punctures. 

Following [8], let us show how Drinfeld's construction in [3] can be adapted for 
X = CP 1 with marked points z a . Consider the space 

Mg\x, Za) = {(parabolic structure on 2 CP \ at z a , 

class of morphisms O CP i — > Ogpi)}, 

the morphisms being considered up to Oc pi -automorphisms; the cokernel of the 
morphism considered in this definition is Ocp 1 - It is natural to consider K^pi Za = 
^CP 1 (^2a=i( Za )) as the canonical bundle in our punctured situation, and then the 
space {classes of morphisms K~ pl — > O CP i}, its mapping n to X N ~ 2 (given by 
the zeroes of a given section) and the diagram 
(1) 

Z 

Po / \qo 

?/ \ s 

M%\CP\z a ) {0 C pi} 

p being the projection on the first factor, and the correspondance Z being de- 
fined to be the set of ((l a ,i),j), l a : line in the fiber of O cpl a ^ z °" ^ m orphism 
C^cp 1 C^cp 1 ' 3- m orphism K CP i Za — > CP i, proportional i,j's being considered 
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equivalent, such that denoting by k : 0^ pl — > Ocp 1 , the cokernel mapping of i, 
there exists a lift j' : Kcpi Za — ■> C^p 1 of j (i.e., we have j = koj), compatible with 
the parabolic structure (i.e., the image of j' at z a should be the line l a ). 

Let us fix weights X a for each a; the P-modules we will consider will be twisted 
by the following line bundles: £(\ a ) on M.q (X, z a ), £>*£( Aq ) on M.^ (X, z a ), and 
PoP*£(\ a ) on -2- For ^ a variety and £ a line bundle on F, we denote (V Y )c = 
£ ® XV ® Let us fix now complex numbers /j, a ,a = 1, • • • , TV, s.t. X^=i A*a = 0, 
E«=i + E«=i 2A a (A a - 1) = 0, £^ =1 f/ a z* + £)^ =1 4A a (A a - l)z a = 0; we 
associate to them the operator on X, 

/ \ n r>o2 \^ \^ 2A a (A a - 1) 

a=l a=l v 7 

and the £> x -module £(a q ),(m«) = V x /V x D (Xa) ^ a) . 

Consider on the other hand on M.q (X, z a ) the twisted "D-module 



N 



, 2 

2''' 



(the conditions on /i a correspond to the relations on the L a , 

N N N .. 

Y,L a =0,J2 L « z * + Yl 2A «( A « ~l) = ef + fe + -h 

a=l a=l a=l 

N N 

a=l a=l 

e = X^=i e ^> e i = ^2a=i z a e ^ a \ analogous relations for /, /i, /i, /ii). We would 
like to show: 

Proposition. — (cf [8].) There is a homomorphism of V-modules 

«*£fa% a )^R(to).P&M 0la) [N], 

which is an isomorphism over it~ 1 ((X — {oo})^ -2 -* — A) (A is the diagonal part of 
(X-{oo})( N ~V). 

(Here we denote, for T a sheaf on a manifold V, by the sheaf {pv)*(F® n ) 
on y( n ) = V n /S n , pv being the projection V n — ► V^ n \) 

Proof. Let us give coordinates to the spaces of diagram (1). 

M%\x,za) ~ G\(G/B) N = GMCPY, 
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choosing the identification CP 1 ~ G/B, 1 1-> ^-1 i ^ ^' ^ ^ ^ i o ) ^' t ' ien ^ 

acts on (CP 1 )^ by homographic transformations. Now A4 B (X, z a ) = B\(G/B) N ; 
after fixing Ocp 1 — > @cpi ^° ^ e the nnes are C(l,t~ 1 ), the t~ x being 

defined up to a global affine transformation. An element of B.om(K~p 1 z ,O CP i) 

is a 1-form u " dz , with « Q = 0. The incidence relation defining Z is 

ti a t a = 0, since the first component of 7' has to be Ua ^_ adz and should 

be regular at 00. The map 7r associates to (w Q ), the solutions (wi) of = 

(counting /c times 00, if this function is ~ c/w 2+k for w — > 00, c 7^ 0). 

Let pi be the natural projection of (G/B) N on G\(G / B) N , then 



PiCM(m«)) = V (G/B)*/ ^2v {G/B)N (L a - fj, a ) +V {g/b) ns1 2 (C). 



Introduce the formal variable 2, then 



E ^Fr 1 + E = «w/w + /w«w + 5 ft w 2 - E 7^° 



Q = l 



2 2 a (2 2 a ) 2 2 2; 2 C 

a=l a=l 



with e(z) = ELi etc., e<«) = £ £ + 2 A a t a , /(«) = = 2(f a £ + 

A Q ). The Radon transform of the P-module generated by the L a — \i a is the "D-module 
generated by the L a — fi a , where 



N f N o\ (\ _ -\\ 1 N 

E ttt 2 + E ,f z J> = *w/w + im*) + - E T^r 

a=l a=l a=l 



e ( z ) = Eli F=~T' anal °g° us formulae for f(z), h(z), = -(u a (-^) 2 + 2(A Q + 
/ (o) = ««" = -2(u a £- a + A a + 1). 
Consider the operator 



N 1 



a=l 
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and let e(wi) = J2a=i w -z e< ' Q ' > ' analogous formulae for f(wi), h(wi). Then 

, tfi , \ w i z a) z 

a=l a=l 



l<af<N Wi - Za Wi ~ Zf3 Wi ~ Z/3 

2 Wi- zp 

1 r g g / 1 

2 8 . 1 . . d . 

+ M Q - ( )( W/3 __) 

- Z a dU a Wi - Zp OU/3 

- £ — {[2(A a + l)/-,^— W 

2 9 / 1 

«aa^(- ^)( A /3 + !)} 



with [a, b] + = ab + ba. Now, 

TV 



tir-(- L -)'v 



E u «(f-) f )up = -2u a -^-( ), 

^ Vdw a / \Wi-zpJ du a \Wi-z a ) 



so the last line gives zero, and the term in [ , ] + gives 

N o iV 



i ^ i i a / i 

2 / . h 2 / v w « a 

, w 8 - 2 a dw a Wi-u a ^ Wi-z a du a \Wi-u a 

a=l a=l 



^ (w, - z a ) 2Ua du a ' 
(the dot denotes the product of differential operators). Then we deduce from 

-r-rTV-l/ \ N 

( f + ^) ELl ^ = ^ S ° d «« = ^(f + ^t) ^ 

^— in- — 



dw; ^ Wj - z a <9-u Q ' 

a=l 
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so the remaining term gives -2 £^ =1 ^-(_l_)( u ^^-) = £^ =1 ( ^_% )2 (^^)- 
Finally, 

L(iUi) = e(wi)f(wi) + f(wi)e(wi) + \h(wi) 2 - — 3 2A Q (A Q 1) ^ 

a=l a=l 

Now, f(w,) = 0, and h(w t ) = -2[^- + with A(w,) = E«.i so 

^) = + ~ it zrrr - it 2XJX ° ~ 1] 



■tWi-Za ^ {Wi - Z a f 

= 1 a=l 



In this way, we have constructed an epimorphism from the (N — 2)-th symmetric 
power of £(A a ),(M a J (restricted to the complement of diagonals) to the D-module 
generated by the L a — p^s (restricted to the complement of the discriminant), and 
so to the P-module generated by the L a — p^s and the action of sfaiC). 

Let us now show that it induces an isomorphism of the sheaves of local analytic 
solutions. Let us start with a local homomorphism of the first sheaf to the (analytic) 
structure sheaf. It is some function (ip(wi))i<i<N -2, such that (d^. + q(wj))ip(wi) = 
0, for all j. We deduce from that relation, using our previous computations, 

(3) Y ( L °~^ =0 , z = l,...,iV-2. 

Wi — z a 
i=i 1 a 

We will consider ip as a distribution on the space of all (tti)i<i<jv, supported on the 
hyperplane ^i=i u i = 0' an d analytic on this hyperplane. We obtain from (3) 

(4) (L a - fj, a )ip = u a (f) + z a u a p, 

<p and p being distributions of the same nature as ifi; indeed, <p an d p can be obtained 
solving a Cramer system (since everywhere on the support of ip, we can find two 
indices a 7^ (3 such that u a up 7^ 0), and L a 's commute with J2 a =i u a- From the 
relation J2a=i = follows that p = 0; from J2a=i L a z a + E«=i 2 ^a(K ~ 1) = 
ef + fe + ^h 2 and EL^ + EL^f^" 1 )^ = 2 (^f + he+ \h\h), follows 

that (Y^a=l u a)^ = (52a=l z aU a )(f>, and (Eli^«a)# = (E! = l4 tt a)f So ; on 

the complement of {(« Q )IEa=i z A = or Ea=i z « u « = or (E! = iZ a ««) 2 - 
(X^2=i w a)(Xl L I=i Z a u a)}, 4> an d will vanish. So these distributions would have 
to be supported on a subvariety of the set of all (u a ) of codimension > 2, which is 
impossible. So we will have 

(5) (L a - fj, a )ip = 0, eij> = 0. 

Since fip = hip = by construction, we have shown that ip can be considered as a 
local homomorphism of the £>-module generated by the L a — p a (restricted to the 
complement of the discriminant) to the structure sheaf. 
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These two morphisms are clearly inverse to each other; if we show that both 
"D-modules have their characteristic varieties supported on the zero section, this will 
prove the proposition. 

For ) 2 (v ) ^ * s c l ear 5 since it is true for £(A a ),(>„)• On the other hand, the 

characteristic variety of R(qo)*PiP* M^)^] is the set of with J2iLi u i = 

and up to equivalence (-u^, £j) ~ (-u^, £i + A), and 



N N c N c2 

^ Z-Zi ~ \^z-z l )\^z-z l 
1=1 1=1 1=1 



This equation gives 



N ^ a2 , \ N 



* - * n r=i(* nr=i(^ - *) s z - Zi n =i(* - *) 



i?, A, £? polynomials. The set {lUj, 1 < z < iV — 2} is the union of the set of zeroes of 
RA 2 , and of oo counted N — 2 — degRA 2 times. Since the w^s are pairwise distinct, 
we have A =const., and degi? = N — 2 or N — 3. Since no Wi coincides with oo, 
deg-R = N — 2; so, B is also constant, and the £j's are all equal; but this is equivalent 
to & = 0. ■ 

2. Separation of variables for the s/ 2 Gaudin-Calogero system 

Let X be the elliptic curve C x /(/ z , with marked points z a q z , a = 1, ...,N. In 
the s/ 2 case, the Gaudin-Calogero system (which plays the role of the Gaudin system 
in the present situation, cf. [4], [9]) takes place in the space 

M.q (X, z a ) = {(£(t,t-i)) parabolic structure at z a given by t a G CP 1 )}/ 

[(t,t a ) ~ {qt,z a t a ),(t,t a ) ~ {t,ut a ),u G C x ],t G C x ] 

= C x x (CP 1 ) 7V /C X x Z N . 

Here £(*!,...,*„) is the bundle on X defined by C x x C n /[(^,^) ~ (gz, diag(tj)£)], for 
ti, ■ ■ ■ ,t n G C x . 

We consider then the space 

M%\x,z a ) = {{£ {t ,t-i),j ■■ £t -> £(t,t-i), par. str. given by t a G CP 1 ),* G C x }/ 
[j ~ Aj, A G C x , (Ma) ~ (qt,z a t a ), (t,t a ) ~ {t,ut a ),u G C x ]; 

it has a natural projection p to .M^(X, z a ). Consider now the diagram 
(7) 

Z 

Pi/ \Qi 
M%\x,z a ) {£ G Pic°(X),u; GPHomf^ 

P / Po\ / Q0 \ 7T 

M%\X,Z a ) Pic°(X) I^' 
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N 

where K x , Za = ^ X (Y1 ( z a)), Po is the projection 

a=l 

class(£( M -i), j : S t -> ^(t,t-i), Par. str.) i-> £ t -i 

(po associates to j its cokernel), go associates £ to (£,u}), tv associates to (£,u;) the 
set of zeroes of u, and Z is the incidence variety, defined by the conditions that u lifts 
to a morphism j' : K Xz ^ <E> £*-i — > £(i,t-i), compatible with the parabolic structure. 
Writing 

_ A e{t- 2 zz~ l ) dz 

the first component of f has to be Y^_i u a t a t[ ZZ -i] — , so that the incidence con- 

dition is Ea=i = °- 

(Fix our conventions for 9- and p-functions: 9(z) = Yli>o(^ ~ Q lz )Yli>o(^ ~ 

q'z' 1 ), p(lnz) = -(!)"(*), so p(r) ~ r" 2 + • • • for r -> 0; we denote /(z) = zg.) 

The lift to C x x (CP 1 )^ of the Gaudin-Calogero operators ([4], [9]) are defined 
as follows: let z be a formal variable, belonging to X — {z a }. We have 

L{z)=e{z)f{z)+f{z)e{z) + 1 -h{zf 

N 9 N 
= L + J2 L^z- 1 ) + 2A «( A « " l)p{^zz- x ) 



N N x 

+ E^ (a) (E /i(a) )(^(^" 1 )) 2 ' 



2 ^ /v,9 

a=l a=l 



where e(z) = E« = i ^-^V ^' = + + E ^ l^^' 

+2 



/(*) = £«=i^£V o) > and eW = ££ + 2^, /<"> = fc(«> 



at 



Let us fix now complex numbers /ii a , a = 0, • • • , AT, with Ea=i »<* = (this 
condition corresponds to the fact that "52^=1 L a belongs to the left ideal generated 
by Yla=i h^)', consider on A4q (X, z a ), the D-module (twisted by the quotient 
£k,(\ a ) of ^- 1 (^fc 2 ) B ®Z=i£\ a , * ■ Ker( S ) - X< 2 ), s : X( 2 ) - X the sum mapping, 
£/c a bundle of degree on X), 

N 

Consider then the operator on X — {z a }, 

d \ 2 N 9 N 

D (\ a ),(p a ) = 2 ( W Q^J ^{WZ- 1 ) + 2 A «( A « - VpQnWZa 1 ), 

a=l a=l 



and the T>x -module 

The fibration 7r has fibers C x ; we will twist inverse images under 7r by the function 

CL a =i( Aa+1 ) ((7 coordinate on the fiber). Also, we will work with k = (multipli- 
cation by 9{t 2 ) k taking us back to this case). We will show that: 

Proposition. — There is a homomorphism of V-modules from the twisted inverse 
image tt*S^^ to R{qi)*p\p* M( Pa )[N], which is an isomorphism over tt~ 
A) (A is the diagonal part ofX^ N) ). 

Proof. Let 

pi : C x x (CP 1 )^ -> C x x (CP 1 ) Ar /C x x Z N 
be the natural projection, then 

a a 

Because of the factor "D C x x (cp 1 ) Ar ^ )> the P-module is constant along 
the fibers of the action of C x . Its Radon transform is the D-module generated by 
L Q - p a and £)a=i where 



AT A N N N _ a 

a=l a=l a=l a=l 

= e(z)/(z) + /(z)e(z) + ifc(z) 2 , 

= E^li M^-2\(!f Za -i\ ^ a K analogous formulae for f(z) and h(z), with 

e (a) = -K(t^-) 2 + 2(A a + l)J-}J {a) = u a , ^ (Q) = "2[u a /- + (A Q + 1)]. 
Consider the operator 

TV a TV a N _ 

L( Wi ) =L -po + Yl ofa'a^La ~ **«) + E (fl (W)) 2 3 ^ £ 
a=l a=l a=l 

Let e(wi) = X^=i H^iW^^T ^^^ e ^ c - Let us compute the difference 



1 ~ r $ 

-L(wi) + (e(wi)f(wi) + f(wi)e(wi) + -h(wi) 2 - A*o - E ^-^(wiz' 1 ) 

a=l 

+ 2A Q (A Q - l)p(lnw^- 1 )); 
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it is equal to 



(8) 



1 N f> f> 

+ 1 E ^z-^^^w^m. 



2 V 

a,/3=l 



The first term of (8) is the sum of (9) and (10), where 



a,/3: 



AT 



and 



(10) - E i ^(W) ( - 2)(A « + 

= f2(A Q + l) 



a=l 
N 



= 2 ( A « + l)b(lnt 2 ) - ^(lnw;^- 1 )] 



a=l 

We have 

AT 



1^ fh, W + 2W„„..-1\^ 



and 

AT 



K e{t*)e{wiZp x y 13 du a { e{t 2 )e{ Wl z^y 
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So 

N 



1 ] ^ a e(t^)e(w i za 1 )e(t-^e(w i za 1 )du a 



N q 

= ^(-2)[p(ln«;^- 1 ) - p(lnt 2 )]u 



du a 

a=l 



The second term of (8) is the sum of (11) and (12), with 

, N v — v . —i\ 9 ,9 , — 1\\ 9 

(11) =2 2^ n(WiZ a ) U oc^—{n\ W i Z l3 )) U P 



and 



a,/3=l 



( 12 ) = ^ E ^(^a 1 )2^^-(^(^ 1 ))2(A^ + l) 

= E( A « + 1 H- 2w; ^^(^ 1 )]} 

a=l 



^2(A a + l)p(ln(«;^- 1 )) 



Q = l 



To compute (11), we express the relation between the and the we have 

SO 

by inspection of the pole at z a , and because of C ff*~ 1 — — ^nr = u a , it follows that 

AT 

du a =u a [—-2^—-{z a w l )], 

i=l 1 



and so 



Wi- — = - > -(z n w, )u f 



dwi ^ 9 a 1 a du a ' 

OL = l 
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Now 



N d 9 d N d 

(11) = 2g-«, i ^-(-(^«,ri)) W/3 _ = 2gp(ln^- 1 W^-. 



Finally 



(9) + (11) = 2p(lnt 2 ) = 2p(lnt 2 )C^, 

a=l 

and 

JV 

(10) + (12) = ^2(A a + l)p(lnt 2 ). 

a=l 

The sum of these terms is 2p(lnt 2 )(C^ + ^2^ =1 (X a + 1)); the term C-^j + 
X^=i(^a + 1) (equal to — ^ is set to zero in the twisted inverse image 

-K*£f^\ , s. It follows that 



AT 



d d J\ 9 

L(Wi) = 2{t 2 TT75 + Wi fT~- + A ( Wi ^ 2 ~ A*0 - X] ^oi^a 1 ) 



a=l 



N 



+ ^2 2 -^a(A Q - l)p(lnty i 2 a 1 ), 



a=l 



since e(tUj) = 0, and = 2[t 2 ^ + fcf (t 2 ) + w^^g- + A(w;)], with A(w;) = 

— X^=i(^a + l)f ( w i z a 1 )- The addition of the term A(tUj) corresponds to the twisting 
by a GL(l)-connection and does not change the PGL(2)-oper, as in [8]. The addition 
of the term t 2 d/d(t 2 ) corresponds to the fact that we are working here with the 
variables (wi,t), which are linked by the relation t 2 YliLi w i = Tia=i Za ( m °d- Q Z )- 
Our statement follows as before. ■ 
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Remark. As in the rational case, one may remark that the conditions on the /i Q 's 
to satisfy the Bethe ansatz equations, that can be found in [7] , can be translated into 
the condition on the projective connection defined by them, to have a single-valued 
solution ip(z) = Y\ O(zai) / n«=i 0(zz~ 1 ) Xa . Thanks to the Leray formulae for Radon 
transformation [1], one could expect the Bethe eigenvectors to be expressed in the 
form 

EN 
a=1 (^+i)^( Wl )...^( Wjv) 

N 

■ 2^(— l) a u a dui A ■ • • A du a A ■ ■ • A duN, 

a=l 

1 < k — J2a=i(^a + 1) < N — 1, the integration being on a suitable cycle in CP N . 
In the general case, this formula should lead to the computation of the monodromy 
of the Gaudin-Calogero system (by deformation of the cycle of integration) . It might 
be interesting to express this monodromy representation directly in terms of the one 
of the projective connection associated to the /U a 's. 
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